Quantum Potential for Diffraction and Exchange Effects 



James Dufty and Sandipan Dutta 
Department of Physics, University of Florida, Gainesville, FL 32611 

Michael Bonitz and Alexei Filinov 

Institut fiir Theoretische Physik und Astrophysik, 
Christian- Albrechts-Universitdt zu Kiel, D-24098 Kiel, Germany 
(Dated: March 17, 2009) 

Abstract 

Semi-classical methods of statistical mechanics can incorporate essential quantum effects by 
using effective quantum potentials. An ideal Fermi gas interacting with an impurity is represented 
by a classical fluid with effective electron-electron and electron-impurity quantum potentials. The 
electron-impurity quantum potential is evaluated at weak coupling, leading to a generalization 
of the Kelbg potential to include both diffraction and degeneracy effects. The electron-electron 
quantum potential for exchange effects only is the same as that discussed earlier by others. 

PACS numbers: 
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I. INTRODUCTION 



The application of classical Newtonian mechanics to materials is generally limited to con- 
ditions of small characteristic quantum wavelengths (e.g., high temperatures, large mass). 
In some cases (e.g., electron-proton systems) there is no simple classical limit due to the sin- 
gular attractive interaction. Still, it is useful to explore possible realizations of an inherently 
quantum description as a semi-classical problem. This cannot be done in general, but exact 
or approximate correspondences can be made for specific properties. The advantage of such 
semi-classical realizations is that powerful classical methods can be employed to address 
the difficult many-body problem (e.g., Monte Carlo integration methods for partition func- 
tions represented in terms of classical actions, molecular dynamics (MD) implementation of 
Newton's equations). 

One approach that has met with significant success is to replace the given interaction 
potential with an effective "quantum potential" in a corresponding classical description. 
The quantum potential incorporates some or all of the important quantum effects in a 
modification of its functional form. In the case of pairwise additive potentials, a quantum 
potential has been defined for equilibrium calculations in terms of the exact two particle 
density matrix for a given pair of particles by equating it to the corresponding classical form 
with an effective potential. In this way, the quantum potential incorporates the quantum 
diffraction effects and other non-classical features such as binding energies. A practical form 
is obtained by a first order expansion of the quantum potential in terms of the given potential, 
leading to the Kelbg potential In the case of the Coulomb interaction, the Kelbg form 
shows a "regularization" of the short range singularity by a smoothing of the potential over 
distances of the order of the thermal de Broglie wavelength. Important applications of these 
potentials include MD simulations for a Hydrogen plasma, and construction of an action for 
the singular Coulomb interactions to allow path integral Monte Carlo (PIMC) evaluation 
of quantum partition functions. More general non-perturbative methods to determine such 
a quantum potential from the two particle density matrix have been explored and tested 
3, 4|. Extensions of these ideas to external forces [sj] and non-equilibrium states also have 
been discussed {g]. 

The most important cases of interest involve electrons under conditions where quantum 
degeneracy can be important. Quantum potentials based on the two particle density matrix 
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do not account for many-electron exchange effects. An important early study of this problem 
was the construction by Lado of a classical system incorporating the exchange effects of an 
ideal quantum gas [7]. The classical gas has pairwise additive quantum potentials chosen 
to give the correct quantum electron - electron pair correlation function. This was accom- 
plished using the second equation of the Born-Green hierarchy Q], solved for the quantum 
potential in terms of the known correlation functions. This idea has been given a more 
practical form with the role of the second Born-Green equation replaced by the hypernetted 
chain (HNC) integral equation approximation 8| relating the correlation function to the 
quantum potential [9|. Subsequetly, the interacting quantum system is represented by an 
extended quantum potential that is the sum of that described above for exchange effects 
plus a regularized real potential of interaction with diffraction effects. The objective here 
is to illustrate the simplest case in which the effects of degeneracy and diffraction appear 
coupled, rather than additive. The system considered is again the ideal Fermi gas, but 
with the addition of an impurity interacting with each particle. The corresponding classical 
system has an electron-electron quantum potential as described by Lado for exchange, and 
an additional electron-impurity interaction with both exchange and diffraction effects. An 
additional Born-Green equation for the electron-impurity quantum potential entails a new 
correlation function for the impurity with both diffraction and exchange effects, as well as 
coupling to the electron-electron quantum potential. This equation is solved for weak cou- 
pling conditions, leading to the Kelbg result in the non-degenerate limit but more generally 
describing coupled exchange and diffraction effects. For the case of Coulomb coupling to the 
impurity, it is shown that the degeneracy effects can be described to good approximation by 
an appropriate scaling of the Kelbg functional form. 

There are many different ways in which attempts have been made to introduce quantum 
effects into classical descriptions, so it is important to clarify the context of the present 
calculations. First, they are among a class of quantum potentials that are based on equi- 
librium properties and pairwise additivity. Their use in molecular dynamics simulations for 
nonequilibrium states and for transport properties are therefore uncontrolled. Three-body 
and many-body quantum effects are not included so the formation of bound pairs may be 
described accurately ^\ but more complex molecular structures are outside the realm of 
accuracy. Representations involving many-body quantum potentials follow directly from 
truncated cluster expansions of the Slater sum and exact field theoretical representations 
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such as a classical polymer action come at the price of considerable additional complexity. 
Quantum potentials not tied to the equilibrium state, such as those from wave-packet molec- 
ular dynamics have a potentially wider domain of applicability, but also entail a new level of 
phenomenology. A more controlled introduction of momentum dependent quantum forces 
from the Wigner representation of the von Neumann equation are specific to each state, 
equililbrium or non-equilibrium, but are still in an early state of exploration. A closely 
related field is that of quantum hydrodynamics. Some of the diversity of issues around 
quantum potentials have been critiqued recently [l^ 

It is a pleasure to dedicate this contribution to Frank Harris - exceptional mentor, col- 
league, and friend to all fortunate enough to have crossed paths with him. 



II. QUANTUM POTENTIALS FOR IMPURITY IN AN IDEAL FERMI GAS 

Consider a system of non-interacting electrons at equilibrium in an impurity field fixed 
(e.g., infinite mass) at the origin. The Hamiltonian operator is 

where Vi^a) is the central potential due to the impurity at the position qc of electron a. 
A caret over a symbol is used to distinguish an operator from its corresponding classical 
variable. The average number density at a distance r from the impurity in the Grand 
Canonical ensemble is 

ne.(r; z, /3) =< n(r) >=\Y. Tre-^^n{r). (2) 

^ N 

Here Tr denotes a trace over a complete set of anti-symmetrized N electron states. Also, 
the partition function Z and number operator n are 

N 

Z{z, P) = Y, z^Tre-^^, n{v) = 5 (r - q„) , (3) 

N a=l 

(3 = l/ksT is the inverse temperature, and z is related to the chemical potential /i by 
z = e^^. Similarly, the pair density for two electrons at distances r and r' from the impurity 
is given by 

neei(r, r'; z, (3) =< (n(r)n(r') - 5 (r - r) n(r)) > . (4) 

4 



Finally, all correlation functions for electron densities at arbitrary positions not referred to 
the location of the impurity become independent of the impurity in the thermodynamic limit 
and therefore are just those for the ideal Fermi gas, e.g. 

nee(|r - r'l ; /?) = neei(r, r'; z, /3) |y=o. (5) 

A corresponding representative classical system is defined by the Hamiltonian 

The " quantum" potentials Vei and Vee are chosen to assure that the classical system preserves 
key properties of the underlying quantum system. A natural choice is the requirement 
that the classical electron density about the impurity nei(r; z,l3) and the classical electron- 
electron pair density nee(r, r'; /?) be the same as those for the quantum system. This 
requires calculation of the classical expressions for rigj and riee for the Hamiltonian 
as functionals of the quantum potentials, equating these expressions to the corresponding 
quantum expressions, and inverting those equalities to find Vei and Vee as functionals of the 
quantum rigj and Hee- Although calculation of the quantum expressions is straightforward 
(but non-trivial for Uei), the corresponding classical calculation confronts the full many-body 
problem due to the pair interactions in ([6]). Lado approached this problem by considering the 
exact Born-Green equations obeyed by the classical forms for nej(r; z, (3) and nee(r, r'; /3) 

Vi^ei (ri) = -(3nei (ri) ViVei (ri) - 1^ j dr2neei (ri, ViVee (r2i) , (7) 



Viriee (^2) = -Priee (tu) V iVee (^12) - (3 J dr^Heee {^1, r2, Ts) VlVee {r^i) . (8) 

These equations are part of an infinite hierarchy, coupling correlations among m particles 
to those for m + 1. For example, ([8]) relates Uee (ru) to the quantum potential Vee (^"12), as 
desired, but also couples it to ngee (ri, r2, rs). In the present context, ngg (^^12) is replaced 
by the known quantum form, but neee r2, ^s) must still be calculated as a functional of 
the quantum potential. Then can be solved for Vee (f^u)- Thus, the difficult many-body 
problem reappears in the need to calculate rieee (ri? 1*2, rs). A similar difficulty is clearly 
present in equation ([7]) for Uei (ri). 

Lado avoided the classical determination of neee by using the corresponding quantum 
correlation function, a much easier ideal gas calculation 7|. Then ([8]) becomes a simple 



linear integral equation that can be solved for Vee numerically. However, this use of the 
quantum expression for rigee introduces a new approximation since ([8]) follows from the 
classical Hamiltonian in terms of the classical form for functional of Vee- There is 

no reason to expect that the classical and quantum forms should be the same. An alternative 
approach [9] has been suggested more recently based on a classical "closure" expressing Ueee 



in terms of Ueeifu) and Vee, the hypernetted chain (HNC) approximation 8|- This is an 
approximation to the classical many-body problem and therefore more self-consistent than 
the Lado approach. In practice, it is found that results obtained by both methods are quite 
close. 

Since ([H]) is determined independently of the impurity it will not be considered further 
here, and Vee will be considered as known for the purposes of solving ([7j). The latter has 
similar problems to that just described, namely determination of the classical form for Ueei- 
In addition, the quantum form for Uei is more difficult, requiring construction from the 
eigenvalues and eigenfunctions for an electron in the presence of the ion. This is similar 
to the problem considered by Kelbg for the two particle density matrix. He simplified 
the problem by considering weak coupling conditions, and the same will be done here in 
the remainder of the manuscript. Weak coupling here means l3V << 1 so that functional 
expansion of Vee, ""-ei ('"i) , and neei (ri,r2) can be exploited. This is described in the next 
subsection. 

A. Weak coupling 

It can be shown from ([7]) that Vei vanishes if = 0, and so can be written 

/3Ve.(r \V) = j dr'G{\r - r'\)(3V{r') + .. (9) 

The dots denote second and higher orders in PV. Similarly, 

He. (r) = He {z, P) + J dv'^-^^ \v=o V{r') + .. (10) 

where Ue {z, (3) is the ideal Fermi gas density. Finally, the classical definition for Ueei (ri, r2) 
for the Hamiltonian ([6]) gives the corresponding expansion 

Ueei (ri, r2) = Uee {ri2) (1 " PVei (n) - (3Vei {r2)) 

- j (ir3/3Vei(r3)neee(ri,r2,r3) + .. (11) 



Substitution of (fTT!) into the second term on the right side of ([7]), and use of ([8]) gives the 
simphfication 

J dY2neei (ri, T2) ViVee (^21) = j (iraVei (ra) Vin^e (ris) + .. (12) 

With these results, ([7]) can be expanded to first order in V giving the desired expression for 
the function Gdr — r'|) that determines Vei in ([9]) to leading order 

~ /^"^ CT?/ X lv=o= neG(ri-r2) + / dr^ (uee (na) - nf) G'(r3-r2) (13) 
dV[r2) J 

where the ideal gas functions rig and riee {r) are 

ue = / dpn{p\ n{p) = {z-\^^'l'- + l) "\ (14) 



(2s+l)(^^ 1 rfpn(p)e^P-^^ . (15) 



Here, s is the spin of the Fermions. 

The response function fi~^bn^i (ri) jbVij^) on the left side of (fT3|) describes the direct 
effects of exchange and diffraction on the electron interacting with the impurity. In addition, 
this couples via the second term on the right to the exchange effects among electrons not 
interacting with the impurity (i.e. a coupling of Ve% to Vee in ([7])). This coupling is essential 
to describe the degeneracy of the background ideal quantum gas. To illustrate this, note 
that for the special case of V{f) constant, Vei(r \ V) -^V since in that case V simply gives 
a shift of the chemical potential. Therefore, in general 



y drG(r) = 1. (16) 



Integrating f|T3|) then gives 



a-i f ^ fSn{r,z,(3 \V) due . f ^ ( ( \ 2\ 

~ ^ j — W{7) — ' ''=° " WkTz " y ' ^""'^ ^""''^ ~ '^''> 

= "e-^^^ jdpn^p). (17) 

The second line follows from the definitions f|T4|) and f|T5l) . confirming that the right side is 
indeed the derivative on the left. Thus, it is seen that the coupling of Vei to Vee is essential 
for consistency with the quantum thermodynamics. 
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It is now straightforward to calculate the response function P~^Snei (ri) /SV{r2) at = 
and to solve f|T3l) for G by Fourier transformation. The corresponding Fourier transformed 
potential Vei{k) from is found to be 

where IlQ{k,uj) is the polarization function for the ideal Fermi gas from finite temperature 
Greens function theory 11| 

/, N , f dp n(p) — n (|p — hkl) . . 

Uo{k,uj) = \im2 / — ^ . , , ^ ,, — , (19) 

^ ' n^o J (27r)3 + + e (p) - e (|p - hk|) ' ^ ^ 

e(p) = |^, n(p) = (z-V^(^) + l)-\ (20) 
and dee {k) represents the effects of coupling to Vee 

dee {k) = - [ rfre^^-^ {uee {vn) - nl) = -^^^^ / dpn{p)n{\p - hk|). (21) 

He J Tleil J 

Note that 1 + dee (k) = See{k) is the ideal Fermi gas static structure factor. The quantum 
potential given by f fTSl) is quite general, applying at weak coupling but for arbitrary degree 
of degeneracy. 



III. COULOMB INTERACTION 



An important special case is the Coulomb potential (e.g. a point ion at the origin), 
V{r) = qe/r, where e is the magnitude of the electron charge and the impurity charge q can 
be negative or positive. In the following the k dependence of the coupling to Vee in (I2T1) will 
be neglected (but not its coupling for k = 0). Then the inverse transform of f|T8|) can be 
performed exactly [l^ to determine Vei{r) with the result 

Ve^ir)-.Vir)Sij,z) (22) 

where A = A/27rh^/5/m is the thermal de Broglie wavelength and the quantum regularization 
effect S{r/\, z) is 

4 r (K^ - ) ^ - K^) ) ) • 

(23) 
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Also n*{x,z) is the dimensionless Fermi function normalized to unity and Si(x) is the sine 
integral 



n*{x,z) = „i ^2 ^ , f / dx \ , J , Si{x)=[ dx'^^^ (24) 
It is easily verified that S{r/X, z) is proportional to r for small r/A 

S(,/A,,)^2.^. AW.^=^^^A_, (25) 

SO that the Coulomb divergence is removed. Also S{r/X,z) 1 for large r so that the 
Coulomb potential is recovered, as required by f|T6l) . Finally in the non-degenerate limit, 
z ^ 0, and the Kelbg result is obtained [l| 

SC-, z) - ^^(^) ^ 1 - e--(^A)^ - 2.'- (erf - l) . (26) 

In the opposite hmit of strong degeneracy, z » 1, (1251) gives 




_,_a^^^||. (27) 



where ji (x) is the spherical Bessel function of order 1. 



A. Representation of degeneracy by scaling 

It is interesting to note that the limiting forms fl26|) and fl27|) are both scaling functions, 
scaled by A in the first case and by A/ Vln z in the second case. To explore the extent 
to which effects of degeneracy can be described by scaling alone, consider the degeneracy 
dependent wavelength A {z) defined in ( !25l) . For small z it approaches A while for large z 
it is proportional to l/Vln z as shown in Figure 1. Hence it is a possible scaling length to 
interpolate between these limits. Accordingly, define S*{r*,z) by 

It follows from (l25ll that this scaling assures that the initial slopes of S*{r*, z) are the same 
for all z . Figure 2 shows the extent to which this scaling captures the effects of degeneracy 
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FIG. 1: X{z) /A as a function of ln{z) (open circles). Also shown is the asymptotic limit proportional 
to l/y^ln{z) (dashed line). 

for a wide range of z. Also shown is the corresponding dimensionless quantum potential 
Vei*{r*,z) = S*{r* , z)/r* . For large and small r* the curves are the same, although there 
are some differences for intermediate values of r*. This is due mainly to the oscillatory 
feature that develops for strong degeneracy (related to Friedel oscillations). However, the 
quantitative effect on the quantum potential in these scaled units is quite small. 
This suggests the approximation for arbitrary degeneracy 

S*{r*,z) c:^ S*{r*,0), (29) 

or correspondingly, the approximate quantum potential 

Ve.(r) ^ V{r)SK{^^). (30) 

Here Sxir / \{z)) is the non-degenerate Kelbg form of but now with A replaced by 
X{z). Thus, approximation (!30|) is a universal function for all degrees of degeneracy. The 
change in length scale with degeneracy can be understood by noting that the characteristic 
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FIG. 2: S*{r*; z) and Vei*{r*; z) as functions of r* for z = 0, 1, 10 and 10^. 

energy defining this scale is not fc^T but rather the average kinetic energy per particle which 
approaches the Fermi energy for large z. 

The above explicit results for Vei(r) are limited to weak coupling. In the case of an 
attractive ion at the origin there are important bound state effects that are not included in 
this weak coupling form. However, it has been shown [4] that such strong coupling effects 
can be included approximately by parameterizing the Kelbg form to fit the exact value 
of n(r = Q,z = 0,/5). The possibility of extending this to 2; > in (|30|1 will be explored 
elsewhere. 

IV. SUMMARY 

One of the simplest quantum systems exhibiting both diffraction and exchange effects 
is a fixed impurity in an ideal Fermi gas of electrons. Here, a classical system has been 
associated with that quantum system by the introduction of two quantum potentials. The 
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first is the well-known pair interaction potential among the classical electrons to represent 
exchange, while the second is a renormalization of the bare impurity-electron interaction. 
The potentials are defined by the requirement that pair correlations for the classical and 
quantum systems should be the same. The classical pair potential is determined entirely by 
the ideal Fermi gas correlation function and describes only exchange effects. The classical 
electron-impurity potential differs from the bare potential of the quantum system by both 
exchange and diffraction effects in a complex mixture of the two. A simple representation at 
weak coupling is given by the familiar Kelbg form for diffraction regularization, but modified 
by a degeneracy dependent length scale. 

Applications of classical molecular dynamics to real systems, such as a hydrogen, require 
a classical representation with quantum potentials representing both quantum effects and 
Coulomb interactions among all particles. Current applications use quantum potentials 
for the electrons that are the sum of an exchange potential Vee as determined here plus 
a regularized Coulomb potential of the Kelbg type for diffraction effects. However, the 
analysis of the impurity problem here suggests that exchange and diffraction are not likely 
to be additive. This is clear from Eq. ([8]) where all information about the quantum effects 
enters via n^e where all effects are mixed (e.g., in the random phase approximation). It is 
only in the sense of perturbation in one or the other that they become additive. 
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